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METHOD FOR THE DETERMINATION OF THE SPANWISE 
LIFT DISTRIBUTION 
By A. Lippisch 



SUMMARY 



There are now several methods "by means of which the 
lift distribution of a wing can he determined; hut when 
attempting to apply these methods in practice, it is found 
that they involve, besides considerable loss of time, ex- 
pert knowledge such as a practical engineer may not always 
he presumed to know. So in order to remedy this fault, a 
new Biethod has been developed which even in first and sec- 
ond approximations yields results adequate for practical 
purposes, * 

As in the majority of the other methods, this method 
is likewise based on the Fourier series for the represen- 
tation of the lift distribution. The lift distribution, 
as well as the angle of attack, is split up in four ele- 
mentary distributions. The insertion of the angle-of- 
attack distribution in the Fourier series for the lift 
distribution gives a compound third series which is of 
particular advantage for the determination of the lift 
distribution. For, employing the Fourier series itself 
with a view of representing the lift distribution, involves 
the calculation of a greater number of coefficients neces- 
sary for compensating the ensuing oscillations. The func- 
tion formulated from the Fourier series and the angle-of- 
attack distribution removes this fault. 3y including the 
angle-of-at tack distribution, the effect of the higher 
terms of the Fourier series on the shape of the lift dis- 
tribution is so small that only a few coefficients afford 
an adequate approach. The method is illustrated in an ex- 
ample and supplemented by a graphical method. Lastly, the 
results of several comparative calculations with other 
methods are reported. 

♦"Verfahren zur Bestimmung der Auf t riebsvert eilung langs 
Spannweite." Luf t f ahrtf or schung , June 17, 1935, p7). 89- 
10 6. 
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I . INTRODUCTION 



The chief aim in airplane development is greater 
flight performance, as evidenced by the efforts'of the 
modern airplane designers to minimize the residual drag as 
much as possible. This "brings the useful wing drag always 
more and more to the fore, Uhereas, in the past, this wing 
drag was, on the whole estimated, it now is common prac- 
tice to define the effect of every single wing element in 
order to make the determination of the best shape of the 
wing for the given design problem possible. Aside from 
the flight performances, there are yet a number of other 
factors vitally affecting the design of the wing shape, 
namely, the desired flying qualities on the one hand, 
through which the pitch, yaw, as well as the roll stabili- 
ty are influenced; while, on the other hand, the designer 
will strive to bring the fundamental principles for the 
aerodynamically best wing shape, into accord with the stat- 
ic requirements. As the lift distribution simultaneously 
supplies the basis for the load distribution in the vari- 
ous operating conditions, both problems may be directly 
connected. Lastly, there are the conditions of manufac- 
ture and service, with a view to simplicity and the spe- 
cial service qualities for a particular purpose. 

All these factors affect the design, and it is the 
problem of the aerodynamic calculation of the wing to de- 
fine these effects of such measures and to reconcile the 
different antagonistic requirements. 

This is being illustrated in several examples. The 
best lift distribution for finite span is, as.is known, 
the elliptical. When other limiting conditions, such as 
optimum moment distribution, for example, are demanded, 
they may be "found from Prandtl*s derivations (reference l). 
It further may be presumed, as known, that a minor diver- 
gence from the pure elliptical lift distribution, as oc- 
curs on the conventional wings, has no substantial effect 
on the induced drag. The effect of change of lift distri- 
bution on the induced drag has been pretty well explored # 

But it is important to establish whether the used 
wing form in extreme positions would reveal vitiating 
flight characteristics as a result of the applied lift 
distribution. Such a study is essential on tapered wings, 
for instance. The tapered wing in figure 1, with constant 
profile, and angle of attack, has dashed lift distribution. 
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Defining therefrom the specific loading of the separate 
wing elements, i.e., the lift coefficient c a , it is 
found that the curve of c a along the span reveals a dis- 
tinct maximum toward the tips. In operating conditions 
with high total lift the "breakdown of . the flow then starts 
in' the vicinity of this maximum. As oven the least unsym- 
metry in the air flow or in the design would induce a one- 
sided "breakdown of the flow at a wing, such a wing would 
"be laterally and direct ionally unstable at stalling. Thus, 
the lift distribution must he modified so as to shift the 

maximum c toward the center of the wing. This is ac- 

a 

complished by twisting the wing tips, and figure 1 shows 
that the intended effect has been actually obtained. The 
result on the lift distribution is small and practically 
zero. as far as the flight performances are concerned. 

As pr eviously ' po int ed out, knowledge of the lift dis- 
tribution under the various operating conditions is an 
absolutely necessary .basis for the stress analysis. It 
was common practice heretofore to employ empirical load 
assumptions and to stipulate that they agree with the ac- 
tual loads. With certain limitation, this empiricism may 
prove true for the straight, nontwisted wing, but for 
twisted wings it leads to erroneous conclusions which may 
end in too small dimensions of the structural parts of the 
wing. The lift distribution for various flight attitudes, 

expressed in lift coeff icient of the total wing c aj _ , _ , 

a total 

is shown in figure 2 for the- wing of figure 1. In the 
various flight attitudes the steady dynamic pressure now 
changes in such a way that a drop in c a is accompanied 
by a marked rise in dynamic pressure. Thus when plotting 
the load-distribution curves -the course of c a t q - it 
is readily seen that the load di st'r ibut ions undergo consid- 
erable changes in the different operat ing conditions. One 
striking fact is that the bending moments in diving and 
inverted flight induce enormous stresses at the tips, with 
the result that spar dimensions based on an empirical load 
assumption, would not be safe. 

It would lead too far afield to touch upon the aileron 
effect and torque. Suffice it to say that the seemingly 
statically favorable wing tip modified according to aero- 
dynamic viewpoints, would lose a large portion of these 
advantages . 

In many cases the smooth contour of the wing is in- 
terrupted by cutaway sections or superstructures. Such . 
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local changes in chord distribution must be reconciled 
with a stipulated lift distribution in- order to avoid sud- 
den lift interferences in the range of normal operating 
conditions. Mutt ray's (reference 2) explanation on this 
subject is that the local chord change must be compensated 
by an opposite angle of attack change, or change of camber. 
The middle piece of the wing with cutaway, in figure 3, 
must obtain higher angles of attack, or profiles with 
greater mean camber. 

The inverse process is applicable to wing-body fil- 
lets or cowlings. The amount of this angle of attack and 
profile change can be implicitly deduced from the lift- 
distribution calculation, and the designer is in position 
to effectuate the desired wing changes. 

The few instances cited, already attest to the neces- 
sity of lift-distribution calculations. When, in spite 
of this fact, the application of this knowledge derived 
from airfoil theory is not common practice with the pro- 
ductively engaged airplane builder, the main reason for 
this lack is attributable perhaps to the enormous paper 
work involved with the conventional methods and in most 
cases also to the lack of personnel adequately trained to 
do such work. 

The practical airplane designer needs a method which 
on first approach affords serviceable results and which, 
without the use of complicated equations, conforms to ac- 
tual practice. 

II. GENERA! PRINCIPLES FROM AIRFOIL THEORY 

Notation 





circulation 


ca i 


lift coefficient 


c v 


rolling- moment coefficient 


cn fl . 


yawing- moment coefficient 


F, 


wing area (a 2 ) 


b, 


span (m) 



N.A.C.A. Technical Homo rando Jo , 778 E 

t , wing chord (m) 

y, abscissa along the span 

v, flow velocity at infinity (ia/s) (flight speed) 
w, induced downwash at wing (n/ s ) 

a, angle of attack referred to original flow direction 
CXi , induced angle of attack at wing 
Be -Tl 

2tttv= - lift increase in 2- dimensional flow 

2 B" 



k 



s 



^4 

_ m j 



p , air density 
q> , angle abscissa of span defined with: 



cos cp = - -j- • 

The lift of the individual wing element is a function 
of r, that is, 

a s p v r 

la the lift per unit length. 
Accordingly , 

r = c » 1 ? 

"because D 

A = c a t v - v 2 

is the lift per unit, length. 

The circulation distribution is a function of the 

span 

T - I(y) 

The rate of downwash at point y p due to Y is: 
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+*/3 ®L 

P 4 n -b/2 y P - y 

Owing to the induced rete of downwash the angle of 
attack at the point of the wing decreases "by a^. It is: 

tan tt4 = ~ 



where, since the angle is small, the tangent is replaced 
hy the angle itself (tan a± = a±) . 

The ahove term substituted for V gives 

9 (c a t) 



a- = - — J ■ — dy 

1 8 "-b/2 y p ' 7 



(1) 



and the lift distribution in place of the circulation dis- 
tribution gives 

The lift coefficient c a itself Is" a function of the 
effective angle of attack. It is 



: a - 2 tt T| a effective^ 



(2) 



= 2TT-n(a-ai) j 

whence the integral equation for the lift distribution 

a - -i- / dy 

8 " -b/2 7 P " 7 



C a t = 2 TT T| t 



(3) 



or 



c a 1 



2 TT T) t 



dy 



+V 2 U"L1L 

4 ~b/2 y P " y 



dy 



♦See footnote, page 7, 



(3a) 
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III. ORDER OF LIFT AITD ANGLE-OF-ATTAC-E DISTRIBUTION 

The lift grading given through c a t may "be visual- 
ized as the sura of elementary distributions, so that each 
elementary grading may "be treated separately and any total 
grading may "be obtained for the desired operating condi- 
tions. These elementary distributions are (see fig. 4): 

1. The normal distribution c a ^normal 1 that is, 
that lift distribution for which angle a is the same at 
any point of the wing; that is, a - constant. This dis- 
tribution is contingent upon the plan form of the wing and 
changes with a, and/or the total lift. 

2. The zero distribution c a t£ , which is produced 
through an equilateral - partly positive, partly negative - 
angle-of-at tack curve* The induced positive and negative 
partial lifts cancel, so that the total lift of the zero 
distribution of a semi-wing is always. equal .to zero. The 
angle-of-attack curve therefore corresponds to an equilat- 
eral' twist. The aspect of the zero distribution is gov- 
erned by plan form and twist, although unaffected by 
changes in total lift, and total angle of attack of the 
wing. 

3. The normal roll distribution c a t^: This quota 

of the lift distribution is due to a rolling motion and 
may be visualized as being due to a straight, twisting 
curve. The pertinent a is directly proportional to the 
distance of the particular wing element from the axis of 

OCn 

roll, that is, ~* = constant. The curve of normal roll 

V 

distribution is dependent on the plan form of the wing 



and varies as the rolling moment, or 




4. The roll-zero distribution c a : This portion 

of the lift distribution arises from a steady roll of the 
wing about its longitudinal axis with corresponding aile- 
ron deflection. The rolling moment due to aileron deflec- 
tion is compensated as a result of the induced rolling mo- 
tion, so that the lift distribution produced by it sets up 

*(lrom p. 6.) The treatment of the problem through fur- 
ther reduction of this integral equation is reserved for a 
future report. 
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no rolling moment. The roll-zero distribution is depend- 
ent on plan form and ailerons, ar..d changes with the aile- 
ron sett ing . 

These four clemontary distributions constitute the 
total distribution of an operating condition. If the 
flight attitude is one without aileron deflection tnc roll 
distribution does not take place; the normal distribution 
and zero distribution in twist make it possible to deter- 
mine- any and all operating conditions within the range of 
straight lift increase. For flight attitudes in stalling 
range",' the data from the straight lift increase are appli- 
cable but only with certain restrictions, although such 
flight attitudes lend themselves also to mathematical 
treatment, as shown ' el s ewhere . If an aileron deflection 
is present it results, for the time, in an additional lift 
distribution which creates the particular aileron rolling 
moment. The amount of this rolling moment governs the 
magnitude of the normal roll distribution. If the roll- 
zero distribution "becomes additive to the particular aile- 
ron deflection, it yields the additional aileron lift dis- 
tribution at the inception of roll. During the subsequent 
roll a moment compensation takes place which renders the 
rolling moment zero in the final attitude; that is, the 
normal roll distribution likewise becomes zero. To this 
the roll-zero distribution, conformably to the aileron de- 
flection is then simply added. 

These four elementary distributions have four equiva- 
lent angle-of-attack proportions. They are: 

1. The mean angle of attack , equal at any point 
of the wing and directly proportional to the total lift. 
If = 0, the total lift = 0 also. 

2. Angle of, twist A, which corresponds to an equi- 
lateral wing twist producing, for instance, positive A 

in the center section and negative A at the wing tips. 
The curve of this angle being the result of constructional 
measures (modified profile or incidence), it cannot do 
changed during flight \inlesr, flaps fitted on the wing can 
be deflected equally (landing flaps, trimming flaps, dif- 
ferential ailerons).* 

*The aileron setting with differential operation corre- 
sponds to the deflection without differential for equal up- 
ward pull of "both ailerons. 
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3. The mean angle of roll cu , which, as previously- 
stated, is directly proportional to the distance from mid 
center, or in other words, 0Lq is always a straight line 

a O 

through the origin, and has the sane value at any 

i/ 

point of the wing. 

4 # Angle of twist in roll Aq , produced as differ- 
ence "between the additional angle of attack due to aileron 
deflection and that mean angle of roll which would create 
equal rolling moment. Consequently, Aq produces an ad- 
ditional lift distribution without rolling moment. 

IV. DEVELOPMENT OF METHOD 
1. General Derivation 



The various lift distributions were computed with the 
well-known Eourier series. Since this distribution changes 
with the scale of the wing as well as the dynamic pressure, 
the semispan is chosen as unit, scale and the dynamic pres- 
sure is put q = i. The conversion .of the obtained results 
to the actual dimensions and speeds is effected through ex- 
pansion with b/2 and q. Thus, 

q - i; I = 1 

and all length dimensions are reduced to unit scale by di- 
viding with b/2. The span coordinate y then, is from 
-1 (left wing tip) over zero (wing center) toward +1 (right 
wing tip). The integral for determining C4 is then taken 
within the : -1 to +1 range. 

Expanded in Fourier, series, the lift distribution is: 

c a t = A x sin <? + A 2 sin 29 + A 3 sin 3<p + ... 

- 2 A n sin n 9 (4) 

wherein the angle of abscissa cp is tied to the span ab- 
scissa through 

cosy - ~ y 

sin cp = yT~- ~y*5""~ 



10 



E.A.C.A. Technical 'Memorandum Uo . 778 



The induced angle of attack is; 

a i = 8~s in^ Sin 9 + 2As Sln 29 + 3As Sin 39 + ^ 

= — S n A n sin n 9 - (5) 

8 sin 9 i 

Then the angle of attack is defined as: 

- [ £ + a-] - SH* [« + !] + 



a 



(6) 



• 4.-U ? - Sil1 ^ 

that is, a variable dependent on the chord distribution 
and the lift efficiency of the individual sections of the 
wing . 

As A n sin nt\ appears in the series for the lift 
distribution as well as in that for the angle-of-at tack 
distribution, the former may equally well "be expressed 
v/ith the latter. Eliminating the last term of the series, 
affords 



. sin £ 
c a t = ^ 

6 + 8 



" + SLJ=_1 m si!LSL m_-„2 sin__2£ 
a + 8 Al sin cp + 8 As sin q> 

, m - 3 . sin 3Cp 1 



_ sin eg. 
m 

8 L 



e + m 



n=m . ™1 



, v di - n A sin n ^ 1 7 \ 

a + 1 J 1 ~8~ A * linlp J C7) 



The subsequent calculations of the lift distribution 
are effected with special forms of this generalized series. 
The advantage in simplified calculation accruing from this 
remodeling rests net only in the fact that for a finite 
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arrangement of terms, only m *- 1 coefficients need to "be 
defined, "but also that the fewness of coefficients in an 
approximate solution itself insures a high degree of accu- 
racy in the approach • 

Ths coefficients A x and A 2 are known through tho 
choice of the operating condition; that is, A x . defines 
the total lift and A 2 , the rolling moment. That is, 

+ i 

C *total F = _{ ° a * dy 
/ c a t dy A x 



Ai = c_ , n JL (8) 
1 a total rr 



+ 1 



+ 1 



/ c a t 7 dy 



-l 



4 « 



wherewith, other than a, the terms of the series consti- 
tuting the major proportion of the "bracketed term, are 
known. Interrupting the series at m = 3 or m = 4 for 
the first approach of the lift distribution, it affords a 

correction for c and the rolling moment, which in 

a total 

turn gives the correction for cc m and cCq, thus neutral- 
izing in part the error introduced through "breaking off 
the series, 

T7hen defining the four elementary distributions, the 
series for the total distribution is split up into four 
parts and each is computed separately. The series for the 
normal distribution contains, for reasons of symmetry, only 
uneven terms; that is, X ± , A 3 , A 5 , A 7 , etc. and so does 

that for a zero distribution. In addition, it is A x = 0, 
because the total lift of a zero distribution was assumed 
equal to zero. The series for the normal roll distribution 
contains only even terms; that is, A 2 , A 4 , A s , etc., be- 
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cause the ordinates of the roll distribution for the left 
and right wings are inversely equivalent. The zero dis- 
tribution in roll contains the even terms A 4 , A s , etc., 
because according to definition the rolling moment and con- 
sequently A 3 . ■= 0. 

The total- distribution is formed as the sum of the 
four elements: 

c a Hotal = c a formal + c a *A + c a *Q + C a l Q A (l0) 

whence the series for the total distribution is the sum of 
the series of the elementary distributions; tiiat is, tne 
coefficients of equal sin ncp terras are additive. Tne 
total lift coefficients are hereinafter written with capi- 
tal letters, namely: A x , A 3 , A 3 .,.A n , those of the nor- 
mal distribution with small letters a x , a 3 , a 5 to a 3n+1 , 
or a a . a 4 , a 6 to a 2n , whereas those of the zero distri- 
bution, and zero distribution in roll are designated with 
a prime, thus: a 3 ' , a 5 '...a , 3n+l and a 4 ' , a s '...a an . 



3. formal Distribution 



The determination of the series is explained in de- 
tail as follows. It is: 



c a * normal 



_s in_cp 

7~, 2m +JL 



n % m (3m + 1) - (3n + l) s in_(_2n_+_l}^ 

Oa + n | 0 g- — a (En+i) S in 9 

i.e., for m=3, for example, 

_sin^_ 

c a * normal " 7 
* + 8 

"a + 4 &1 + 4 3 "sin 9 4 ag sin «? _ 

Since the total lift in normal distribution and con- 
sequently . a a .is merely a scale factor, it is: 
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5 a normal _ 



sin 9 
6 + ff „ 



a u m ra - 1 sin 3 9 
ai" + 4 + 4 a x "Tin~cp^ 



m - 2 a 5 sin 5 cp 

H* : 7 XT"* - " H" • • • 

4 a x sin 9 



(11) 



We first approximate a m /a 1 , because 



/ 1 c a * normal 



TT 



Integration of the right side of the above equation 
gives as first integral: 



a. 



+ i 



m nil 



sinjp_ 

t 2m + 1 

e + — 8 — 



ay 



The next integral is: 



+ i 

m - 1 a _ . 

J- = — T— / 



4 



sin 3 <p 

a i -i ? ■ 2mJh_l 
b 8 



ay 



and the other integrals: 

t _ ra - n a aa+i + / ii*_i!. n dv 

211+1 4 a i _I | + ?n_+_l y 

Bearing in nind that i + differs only slight- 

ly from a nean value and that, besides, ^JlJH a s n+ 1 is 

4 ai 

snail, it is found that the integrals from J 3 to J En+ i 

are alr.ost equal to zero, "because 
4-1 TT 

/ sin (2n + l) Cp dy = / sin (2n + l) 9 sin 9 dcp = 0. 
Accordingly, a close approach of a n /a 1 is: 

5a - 2 n , v 

_i £ + S.e-t-1 7 
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The integral is 'best evaluated "by graphical method, 
the numerical calculation being in most cases quite diffi- 
cult. Computing this integral for the first three steps 
m = 1, 2, and 3, the cc m /a 1 values are found to vary only 

a few percent, so that the mean value of these three steps 
is sufficiently accurate- 
Even a mean value for t itself affords a satisfactory 
first approach. Then, 

& x s mean T 8 

From the determination of cuj a. x follows the first 
approximation of the normal distribution, "by "breaking off 
the scries expansion at the third term. Then, 



c a f ormal ®_ = JJ__Zjl 

t + | 



On 1 

— + A 

a, 4 



(13) 



This first approach already affords information a "bout 
the characteristic "behavior of the normal distribution and 
is perfectly satisfactory for design calculations, the er- 
rors for smooth wing contours ranging between 3 and 5 per- 
cent . 

To determine the higher coefficients, we merely com- 
pute that portion of the lift distribution formed as dif- 
ference between the normal distribution and a pure ellip- 
tical distribution for equal total lift: Forcing a pure 
elliptical distribution on the wing by changing the angle 
of attack over the span, that is, writing 



c a t = a x sin cp = a x J 1 - y 2 
the angle- of-attack distribution is given through: 



Elliptic ~ ^ V" ' 8/ 

The difference between the constant oc m and the el- 
liptical angle of attack represents the twist pertaining 
to the difference in lift distribution between normal and 
pure elliptical distribution. The differential grading is 
a zero distribution which changes with the total lift. It 
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is called elliptical- zero distribution. It is: 

c a t normal = c a t elliptic + c a * A elliptic 
1 ikewi se 



k) + A 



ellipt ic 



To define the elliptical-zero distribution from the 
elliptical twist, we form: 



A _ellip_t ic. ?m „ 



c a t A 



e + 



(14) 



a - -elliptic = s 



\ + 



? ^e lliptic + is sin 3 g 

2m 4- 1 ! aj, 4 a x sin cp 

— g— L 



m - 2 ag_ sin 5 cp 
+ — $ £3. TIFT" ' " * 



4. a, 

In first approach, we have: 



la. sin (2m - l) cp 

g cr* i »_ r , 



sin cp 



(15) 



a elliptic 



v 2 



^elliptic 



(16) 



How, to what extent does this approach conform to ac- 
tuality? Subtracting % + 1 from the series for a m , leaves 



^ell i £ t ic = H l^Ll 



a. 



sin 



3 

1 + ? 



a c sin 5 cp 
+ ,1 



+ *S2±± ill±l*l±l± 
a x sin cp 



a x sin cp 



2m + 1 

f + —j 



e + 



+ . . . 



which reduces to 
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"elliptic 
a i 



sin 9 



■J e + 



a. 



sin' 3-'q> + 



a , 



sin 5 Cp 



+ ~- sin 7 9 
a l 



+ —Sti- sin (2m + l) <p 



f + ?.m_t _1 
5 + 8 



The first approach differs from the true series for 
elliptical zero distribution only in the added quotient 



t + 2m_+ 1 
& 8 



i + 



8 J 



Therefore the first coefficients 



etc., may be closely approached from 



Ellip tic 3ln ; 



when 



defining the first Fourier series terms of this approach 
with the aid of harmonic analysis. That is, we write: 



a x 



2m_+_l 
8 



n - c 3 sin 3 9+c 5 sin 5<p + •••• +c n sin n<P 

e + | 

e + 

The quotients 

a mean value with respect to coordinate cp, the desired 
coefficients of the series may "be approximated at 



differing hut little from 



^3 



-= c 3 



l 5 ~ 



a. 



= c 5 



c + - 

5 + ej 



mean 
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A higher approach for 5&JL£&iIiliM£. is then 
*• * a i 



ob- 



tained when inserting the thus-established coefficients 
into (15) such as: 



^a_ t ._^ellip.tic 



r>in 9 



+ I 



^llipt ic + 1 a a a in_3_? 
a a 2 a x sin 9 

1 a s sin 5 CP 



+ 4 a x sin 9 _ 



. .(17) 



To ascertain the decree of accuracy of this approxi- 
mation, we effect various reductions, so that: 



sin 9 



-Qll ial ic + 1 a 3 B in 2 9 1 sin 5 9 
2 a x ""sin 9 4 a x sin 9 



a • 



= a 3L 
a. 



sin 3 9 + 



"5. 



a 



in 5 9 + - 7 - sin 7 9 + - s 



+ sin 11 cp 

a i 



+ . . . 



+ »|s±i. sin (2n + i) 9 



f + 2m_+_l 
& 8 



e + 



8 



Again evaluating this approximation by means of har- 
monic analysis, yields the improved Fourier series coeffi- 
cients d 3 , d 5 , d 7 , etc. After forming the mean values 

I + ?jn_±_l 

for =• the calculation gives the improved approach 

fc 8 

of the desired coefficients of the elliptical zero distri- 
bution at: 
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&3 




d 5 



±1 




n 



'-»mean 




J mean 



This method may "be continued until the desired degree of 
accuracy is obtained. 

It is seen that for engineering purposes in which accu- 
racy need not exceed that of other aerodynamic principles, 
the coefficients may he interrupted with a 7 /a x . In many 
cases even lower approaches are sufficient, when the wing 
contour differs only slightly from the elliptical. 



That is, the distribution produced "by equilateral 
twist. The angle-of-at t ack curve is given and must he 
separated into mean angle of attack and pure twist. In 
other words, the true angle of twist must he referred to 
the direction given through the zero-lift axis of the 
whole wing in place of the arbitrarily assumed body fixed 
reference plane. This separation is effected by bearing 
in mind that the desired zero distribution must comply 
with 



3, Zero Distribution 



i 



/ c a t A dy = 0 



o 



The first approach is: 



LT.a.Ca. Technical Monorandum So. 778 



19 



(18) 



In view of A -f £ "being given rather than A, 
whereby c is the angle of zero lift direction and refer- 
ence plane, we determine € from the integral 



+ i 



/ (A + e) d € / , d 

t + i o e + 3 



8 



8 



From the first approximation of the desired zero dis- 
trict ion, we subsequently determine in similar fashion as 
for the elliptical zero distribution, the first coefficients 
of the Fourier series through harmonic analysis. 

With the zero distribut ion as: 



-a 



; A - a 3* sin z 9 + a s' sin 5 <P + ... 



it is: 

A := a 3 1 

3 . sin H 1 



e + 



+ a s ' 



s i n 5 <p 
"sin~<P 



+ a i §.in_.l2m_+_ll_2 
sm+i sin 9 



f 4 2m_+ 1 
" 8 "J 



so that the first approach of the zero distribution corre- 
sponds to the term: 



A 



sinJP 



sin 3 cp + a 3 ' sin 5 <p 



11 

i + 



+ ... 



~p + 2m_+_l 
8 

~ 1 + 1 

& 8 



With the use of the mean values of 



e + 



2m_+_l 
8 
~3 
8 



t+i 



v;e 
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approximate the first coefficients and then form a higher 
approximation which may "be improved as desired through 
further analysis. The calculation process "being in prin- 
ciple the same as for the elliptical zero distribution, it 
needs no further explanation. In the closer approxima- 
tions, G must eventually "be improved through control of 
the zero lift . 



4. Normal Roll Distribution 



Normal roll distribution produced, according to our 
definition, through straight, contrary wing twist; with 



6m 



as angle of roll, it is: 



cc, 



(19) 



The normal roll di st r ifcut ion is expressed in Fourier 
series with even terms: 



c a = a 2 sin 2 q> + a 4 sin 4 Cp 



+ a 6 sin 6 9+ ... + a 2m 



sin (2m) <p 



The angle-of-attack curve is: 



sin 2 <p f 6 , 2l , sin 4 <P V t , 4[ , 



+ a 



2m 



sin_(_2m^9_ 
sin <P 



i + 



L. 



2ml 
"8"J ' 



whence : 



'a 



, _ sin cp r , m - 1 sin 2 cp 



4 a4 sin cp 



• m_-_n sin (2n)jM 

+ 4 a ^n sin 9 " J • • • • 



(20) 



As the normal roll distribution changes with a 2 , we 
^a ^ 0 

determine — — giving for the rolling moment: 

a 2 
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+ 1 



Integration of the rolling moment from the above term 
for c a tq gives: 



a 



m - 1 

2~~ 



.3 /J _ „2 



+ 1 

. y" V 1 - y* , 

' TTT ay 

-i 1 + 4 



The next integral is: 



+ i 



t _ n - 2 a* > sill 4 <p 



and the other integrals are: 



sin (2n) Cp 



*y 



Then, / y sin (2n) 9 dcp = 0 for n > 1. 

- i 

Moreover, as % + | differs very little from the 
mean value, the integrals J 2 to J n are almost equal to 
zero. Thus, a serviceable approximation for the angle of 
roll obtains from: 



^ m _ 



U 

_8 in j- 1 

, yS i- «L- dy 

o 5 4 



(21) 



The integral of the denominator is defined graphical- 

ly. The mean value of — — obtained for different m 
is within mathematical accuracy* 

In first approximation the integral is determined by 
introducing a mean value for J, Then 
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2 



ft + y 



Having determined — ~~- the first approxinat ion of 



the normal roll distribution obtains at: 



L- a 3 



— "?l TTT" 



(22) 



The other coefficients must be determined if the 
higher approaches are derivated. \7e again separate the 
elliptical roll distribution, given through 

^J-Sellrptic S|h39 



from the rest, a zero distribution in roll, which changes 



as a 



It is: 



(c a tn) A C a tn 

— « — S. = _~ — * - sin 2 cp 

a 3 a 2 

The angle-of-attack distribution pertinent to 
(c a t q )A, is: 



(23) 



a 2 a 2 " sin 9 

A Q /y 



; i 

* + 4 



or, when using -A — for forming the other approxinat ions 



a 2 a 3 



. (24) 



On the other hand, it must: 



-q/-L = 2 a * sin 4 cp 
a 3 



a 3 sin 2 <p 



a s sin 6 
" 2 a~ sTn~2" 



- 2 



Ian ll n _L 2 
a-, sin 



I H + 4"J - ' • ' 
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The first approximation of the zero distribution in 
roll being defined as 



(c a tq)L 



a 2 j " e +" f 



(25) 



the coefficients are computed from: 



i-v J- 



- sin 4 cp 

a 3 



+ gj- sin 5 cp 



t + 2 



"f" • • t * 



+ sin (2n) cp 



£ 4. .ni 
_ _ 4 

i 



e + 



by again averaging 



£ + f 



and then correcting the coef- 



ficients obtained from harmonic analysis. The higher ap- 
proach for the desired elliptical zero distribution in 
roll, then, is: 



a s 



'TT-. ys 

XT}— 



2 a P sin 2 ? 



(m ~ n) a gn sin (2n) 



sin 2 cp 



(26) 



and the coefficients in 1 approach are determined 
through: 
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e + J 



A^/y ( m - 2) a__ sin 4 9 



L a 3 



_4 

a 3 sin 2 9 



(n - n) a__ sin (2n)9 
_ __ ________ 



= *_ sin (2n) 9 + 2 _J*& sin (2n) 9 

n=s a 3 n=m+i t ' 2 



e + I 



To illustrate, for tn = 4, it is: 

y y i - 'y ; 



. 4 



A^/y 2 a 4 sin 4 9 1 a ? sin 6 9~ 



2 a s sin 2 9 2 a 3 sin 2 9_ 



= ~- sin 4 9 + _!*■ sin 6 9 + -J 4 sin 8 9 



01 n 

+ ----- sin 10 9 



t- ---- sin 12 9 
a 3 



+ -ii sin 14 9 
a 2 



e + 



t 1 4 

1 * *J 



The process is continued until the desired accuracy is ob- 
tained # 



5. Zero Distribution in Roll 



The roll distribution due to aileron deflection is 
computed on the basis of a distribution composed of a nor- 
mal roll and zero distribution in roll* As a result, the 
additional angle-of -at tack curve due to aileron deflection 
must be separated into the normal proportion and a twist 
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proportion (fig. 5). As tlie additive zero angle of attack 
due to aileron deflection is proportional to the angle of 
aileron deflection, we plot the ratio of both angles, that 
i s : 



a 



whereas , 



or 



€ = 
G 



cc, 



n 



+ A 3- 

7 



(27) 
(27a) 



quo t a 



Thus tlie problem reduces to finding the constant 

r." „ . e 

for a given - , 

y 



a 4m 



y 



This separation is effected 



on the "basis that the rolling moment due to C equals the 



rolling moment due to <x 



^m 



In first approach, we have: 



_ 



(28) 



with a rolling moment: 
+ i +1 



/ ( 



-1 



t) y dy = / £ y -^-- y 



— e— dy 

+ i 

!a_ / y 

2_i 



ST"! - rrS 



t + 



dy 



Contrariwise, as 

+ i 



/ (c a *)o y dy = - a 3 j 



a< 



is defined with 
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Having, for computing the normal roll distribution, 
defined' the mean angle of roll referred to a 3 , it af- 
fords a Q n /y' 50 ^at the twist follows from the differen- 
tiation of angle roll Aq/jt. 

Formed as "before, the first approach of the desired 
zero distribution in roll then is: 



[(c a t) 0 3 = v "7 T (o0) 



from which the higher coefficients are obtained by harmon- 
ic analysis; a 2 being determined as first approach, the 
zero rolling moment condition must be controlled for higher 
approaches and further corrections effected on Ctq / y and 

V 7 ' 

V, GRAPHICAL ANALOGY 



The mathematical lift distribution may also be re- 
placed by a graphical method having the advantage of clear- 
ness as well as lucidity of the mathematical treatment of 
the problem. The . princ iple is, briefly, as follows: 

Plotting c a t/a against angle of attack a/a x af- 
fords for such lift distributions as formed from the first 
two terms a 1 /a 1 and a 3 /a x the relationship portrayed 

in figure 6. The c a t/a x values pertaining to certain 

points on the span lie on straight lines passing through 
CU 1 Co t 

point — = - ~r> — — = 0. For in view of 
M a ^ 4 a i 



= sin <p + ^ sin 3cp 



3 



c„ t 

it is for ~- — = 0 , 

a i ... 



a_ = _ I 

a 1 ~ 4' 
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With this a desired lift distribution can "be graphic- 
ally obtained in first approach, as follows: 



Find, for a given chord distribution and angle of at- 
tack, the lift distribution giving in first approach a 
certain total lift. In the spatial system of coordinates 



a x " 



(fig. 7), we plot the lift distribution 
conformably to p\ire elliptical distribution 



V ai 



a i b / 



Then we draw the straight 



and establish an 



a x " ~ 4 

area -through the elliptical lift distribution and this 
straight line, in which all lift distributions are con- 
tained which give the total lift 

5a .* 

ai 



dy = + 



rr. 



and can bo formed through the first two terms of the Fou- 
rier series. Fresuming that the exactly possible lift 
distribution in first approach can be interpolated between 
the exact lift distributions, the problem consists in 
finding that curve on the area whose projection marks off 
the given angle of attack from the a/a x to y plane, and 

whose content in c a t/a x to y plane is rr. This construc- 
tion is illustrated in figxire 7 for a normal distribution 
with constant angle of attack* It is readily proved that 
this first approach of the lif t dis t ribut ion is conforma- 
bly to the derived first approach of the lift distribution: 



a 



jr: 



Closer approaches can be determined in like manner in 
the 3-coordinate system by defining, say, coefficient 
a 3 /a x for the second approach (harmonic analysis, first 
approach). Then, however, the track for the area of coor- 
dinated lift distributions assumes the function: 



a 

aT 



track 



as illustrated in figure 8, 
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Further approaches may be formed in analogous manner 
after the corresponding coefficients have "been obtained 
from analysis. Obviously, such a continuation of graph- 
ical determination is of no practical significance he- 
cause, in most cases, the plotting accuracy is insuffi- 
cient to present the still possible improvements of the 
solution. The method is, above all, suitable for intelli- 
gibly showing the method and for effectuating graphically 
the determination of the first approach. 

There is still another graphical method which is even 
more simple. (See fig. 9.) On three planes intersecting 
in a straight line are the functions f (y) , 3P(y), and 
°P(y) together with the stipulation that the points of the 
three functions pertaining to equal y values, lie on a 
straight line. Then, we have, according to figure 9: 

sin (a_+ (0 f(y) 9(y) 

On the other hand, forming the relation for the first ap- 
proach of the lift distribution affords: 



which for 



_2tt v, t /l _£_y 2 _ 

yr - - y 2 + 2tt t) t 



3 



2tt r\ t = f(y) 



y 3 = q»(y) 



cosjt _ 3 
cos P 8 

sin (a + p) _ a_ 
cos P a x 



1 
4 



gives 



This form affords higher approaches when inserting 
the relations between length ordinate y and angles a 
and P as new variables. This method can be readily de- 
veloped so that a mechanical drawing instrument can be 
used for defining the lift distribution. 
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The above-cited derivation for tho graphical determi- 
nation of the lift distribution applies only to normal and 
zero distribution, although there is a similar method for 
the roll distribution. As hereby the rolling moment is 
held cohstant rather than the roll distribution itself, 
the rolling-moment distribution is determined through pro- 
jection. Figure 10 shows the construction for defining 
in first approach a normal roll distribution, where it is 
readily seen that: 



The track is here also a straight line through -1 — = 



+ Farther approaches are deduced in the same manner as 

above, by defining the higher track functions. 

For a wing whose chord or angl e-of-at t ack curve re- 
veals corners so that dt/dy or doc/dy becomes unsteady 
at tofT? points, the method repeated here yields corners in 
the lift distribution in the lower approaches. Since, in 
this case d(c a t)/dy also becomes unsteady in these points, 
the integral for a± gives the absolute term co ; that is, 
the approximation at such points gives a mathematically 
wrong result. Even so, this has no effect on the accuracy 
of the final result. 

To avoid such discontinuity points, the corners in 
chord or angl e-of-at tack distribution can be rounded off, 
as is, in fact, common practice in airplane design. 



The wing is a normal wing of conventional design as 
seen from the semiring illustrated in figure 11. The per- 
tinent data are appended in table I. The efficiency of the, 
lift increase is r\ = 0.89, corresponding to a 1 2 to 14 



graphically evaluated for m = 1, 2, and 3 from the fol- 
lowing table, after which <X m /ai was defined according to 

(12). The compilation appended in the table reveals the 
discrepancies so small that the mean angle of attack can 
be determined within an accuracy of less than 1 percent. 



a 2 L a 2 2 J | + 1 




VI . EXAMPLE 



percent chord thickness* Integral / 

o 




dy was 



f 
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m. 



1 

2 
3 



TABLE II 



r 1 _ _____ 

* P , 2m + 1 
o i + _______ 



it a y 



0.8272 
. 6543 
. 5420 



_m 

ax 



__ 
ax 



0.6995 
.7000 
. 6990 



= .6995 



J mean 



th 



With 0. 695, the value of C^/aj. defined as the 0 
approach over the mean value of i , is also sufficiently 
accurate for rough calculations. A control check for the 
lift increase of the total wing must , of course, give a 
flatter rise than that of an elliptical wing with equal A 
and T). The following tabulation gives the established 
values . 



3c £ 

3a 



ell 



2n ri T---0- - 4.595 

1 A + 2r\ 



TT 



3a 



_ 4.588 



Jeff 



2 T & 



The plan form of the wing "being not much unlike that 
of an elliptical contour, the reduction in lift increase 
is negligible. 

Following the determination of the elliptical ~ 

ai ell 

the four sta.ges of approach of the elliptical zero dis- 
tribution were calculated. The figures are included in 
table I. The determination of the coefficients from the 
harmonic analysis of the individual approaches is readily 
seen in table III in comparison v/ith the coefficients ob- 
tained from the individual stages. 



TABLE III 



m 




_ a _.__„.. 






»i 


a_ 


1 


+0.00 69 


+0.0312 


+0.0022 


2 


+.0063 


+.0314 


+.0020 


3 


+ .0062 


+.0312 


+.0025 
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The existing discrepancies are within mathematical 
accuracy and the coefficients obtained frora averaging 

8 — f e t c , themselves prove sufficiently exact. The 

course of the essential functions is shown in figure 12. 

In order to "bring out the accuracy of the calculation 
for determining a conventional flight altitude, the four 
approaches of the normal dist ributions were reduced to a 
flight condition of c a to tal = 1 *° and tlie lift coeffi- 
cient defined therefrom. The figures are appended in ta- 
"ble IV and illustrated in figure 13. It is clearly seen 
that even the first approach reproduces the characteristic 
form of the lift distribution quite implicitly. The im- 
provements through stages of approach are so minute that 
graphically it is "barely possible to give the fourth ap- 
proach, and in the c a coefficient itself the discrepan- 
cies scarcely exceed the errors involved in the measure- 
ment of those coefficients. 

The calculation of a zero distribution is omitted "be- 
cause it is practically identical with the elliptical zero 
distribution, and we proceed to the roll distribution due 
to aileron deflection. To this end we first establish the 
normal roll distribution of the wing. The figures are 
compiled in table V. It is seen that the values in third 
approach are already sufficiently convergent. For simplic- 
ity's sake °^JJL was merely computed for m = 2, accord- 

ing to (21) at -1.6412 (fig. 14). 

This brings us to the aileron effect itself. The 
shape and size of the aileron is seen from figure 11; the 
calculation process is given in table VI. The aileron ef- 
fect for varying percent of aileron chord was determined 
according to the curve shown in figure 15. The curve was 
defined from experimental data. The reading gave £ = 

-g-p-- for different points of the aileron, while the form- 
ing of the integral: 



> 7 i + i 
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gave the value of a 2 which renders possible the separa- 
tion of the additive aileron setting angle into normal an- 
gle of roll and zero angle in roll Aq/y. From the latter 

the pertinent zero distribution in roll Was then obtained 
in progressive approaches, They are 3ncl\ided in table VI, 
Figure 16 illustrates the course of the essential func- 
tions. As anticipated, the convergence is not as appreci- 
able, due chiefly to the unsteady jump in angle of attack 
at y = 0.40, For this reason it is important to so com- 
bine the mathematical result that the differences in the 
approximate determination of the aileron lift distribution 
are presented with reference to an ordinary aileron set- 
ting angle. Choosing p =+10° affords the figures for 
the roll distribution with 10° aileron setting, given for 
the first and fourth approach in table VII, 



TABLE VII 









c a 




7 


® 


_®__ 


® 


© 


0.00 


±0.0000 


±0 .0000 


±0.000 


±0.000 


.10 


• +.0037 


+.00 53 


+ .012 


+ .017 


.30 


■■ +.0075 


+.0115 


+ .024 


+ .0 35 


. 30 


+.0107 


+.0183 


+ .036 


+ .0 53 


.40 


f +.0134 
\ +.0 834 


+.0269 
+ .0685 


+ .049 
+ .303 


+ .098 
+ . 2 40 


.50 


+.. 0 8 48 


+.07 63 


+ .336 


+ .302 


.60 


+.0838 


+.0813 


+ .364 


+ .353 


.70 


+.0818 


+.08 38 


+ .394 


+ . 40 4 


.80 


+,0771 


+ .0786- 


+ .417 


+ .42 5 


.90 


+ .0G57 


+.0675 


+ .417 


+ .422 


.95 


+.0549 


+.0543 


+ . 445 


+ .440 


.975 


+.0 444 


+.0435 


+ .488 


+ . 4 57 


1.000 


±.0000 


±.0000 - 


i 

1 





The corresponding lift coefficients were determined 
therefrom. The course is seen from figure 17, The errors 
in lift coefficient are minute. To remain within the lim- 
its of accuracy of the measiired aerodynamic coefficients, 
the second approach for computing the aileron effect would 
have siifficed. 

In conclusion, it is pointed, out that an arithmetical 
scheme for a complete calculation of lift distribution 
conformable to the depicted method is being published else- 
where . 
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VII, COMPARATIVE CALCULATION 



In order to prove the practicability of the analyzed 
method, we made some comparisons with other known methods, 
The normal and zero distribution of a wing, computed by 
L.otz's method, was published by S. Hueber (reference 3). 
This is the wing of the "Flavag III" performance glider, 
the computed data on lift distribution having been kindly 
placed at our disposal by the A.V.A. ( Gottingen. Figure 
18 compares the mathematical results for a normal distri- 
bution of the wing| figure 19. the effect of twist. (See 
figs. 5 and 6 of reference 3.) Admittedly, in this case, 
it is not a pure zero distribution because the total lift 
is negative. But, to make possible the use of the numeri- 
cal values , a conversion was forgone* The comparison re- 
veals more pronounced discrepancies between Lotz's figures 
and our own,* but only negligible differences for the 
twisting effect. 

In order to explain the difference between the two 
calculations, we defined the induced Ct^ for both lift 
distributions with our own method* The results are tabu- 
lated in 



TABLE VIII 



y 


a i 
x Lotz 


a i 

Lippisch 


0.00 


+0.0301 


+0.0316 


. 30 


+.0355 


+.0361 


. 50 


+.0404 


+.0378 


.70 


+.0395 


+,0388 



The determination of cc^ involves a graphical inte- 
gration of the curves obtained from the calculation rather 
than an analytical method* 

Next, in order to ascertain the quality of the approach 
of the pertinent lift distributions, the local angle of at- 



*This example was carried through according to the model 
calculation of the report: Aerodynamic Calculation of Air- 
plane TTings (DPS). This model corresponds to the second 
approach of our method. 
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tack is determined at 



a .-n ._S + ai 

2tf ri t 



from c a t , ai and the given chord t. 

Inasmuch as the normal distribution proceeds from an 
angle of attack constant for all points of the span, the 
computed a values must agree if the solution is correct • 

Table IX shows the results for "both methods, while 
figure 20 illustrates the inversely determined a distri- 
bution for both methods • 





TABLE IX 




y 


a Lotz 


a Lippisch 


0,00 


+0 .2499 


+0.2 559 


.30 


+.252 6 


+.2 564 


. 50 


+ .2 300 


+.2565 


.70 


+.2 593 


+.2546 



This ta"ble discloses that the discrepanc ies are due 
to the fact that Lota's method resulted in more inaccurate 
figures. Another method of comparison would he to deter- 
mine the chord distribution from the mathematically de- 
fined mean angle of attack. The rean angle of attach of 
Hueber's calculation being unknown, the comparison was ef- 
fected on the basis of the mean angle of attack of our own 
calculat ion . 



0.00 
.30 
.50 
.70 



actual 

0.2333 
.2 333 
.2185 
.1836 



TABLE X 
t 



Lot z 

0.22 69 
.2298 
,2225 
.1917 



T 



Lippisch 



.2338 
.2187 
,1877 



The results are collected in table X and figure 21. 



Summing up, it is seen that the lift distribution de- 
fined from the second approach of our method already yields 
very practical results, whereby the determination of the 
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second approach - mathematically - can "be carried through 
very quickly. 

Figure 22 illustrates the results of a comparison 
with &lauert f s method for the zero distribution of a ta- 
pered wing with twist, once, according to Glauert's four- 
point method, and then according to the model for deter- 
mining the second approach of the previously cited publi- 
cation. The accord is within accuracy of calculation. 

The check of the calculation method for determining 
the aileron effect was made on an example cited in I. 
Lotz 1 s report, Part III (reference 4), in which the coef- 
ficients of rolling moment, induced yawing moment, etc* , 
are determined for a rectangular wing with variable aile- 
ron length. According to Lotz, the coefficient of rolling 
moment £, is defined as 



TT b 



a c * Fq a q 



wherein 




a 



m 



and 




J b— >co 



= TT 



Lotz T s coefficient of induced yawing moment i is: 




with 



3- 



(a ff = total angle of attack.) 



We determined the first approach of the aileron-lift 
distribution for Lotz ! s constants at 
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and the nornal di strict ion additive' for defining the yaw- 
ing moment from the first approach at 



The calculation of the distributions conformably to 
these relations having been explained elsewhere, requires 
no further elucidation, 

• Our coefficients were obtained through graphical inte 
gration of the particular function. The results of the 
comparison are shown in figures 23 and 24 for the roiling 
and yawing moments of the wing with 1:5 aspect ratio. 

These graphs were taken from Lotz ! s report (figs. 2 
and 3) and our computed points added. Even the first ap- 
proach is already seen to be in close agreement with Lots 1 
curve, while the dashed curve of the determination accord- 
ing to Wieselsberger-C-lauert , discloses greater discrepan- 
cies. 

The local errors of the first approach of the method 
elucidated here, have no effect on the resu.lt s of calcula- 
tions analyzing the effect of the wing as a whole. So 
that even the investigation of aileron effect of varying 
contours and aileron forms should afford practical results 
with the simple relations of the first approach of our 
method. 




Translation by J. Vanier, 
National Advisory Committee 
for Aeronaut ic s - 
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Figs. 1,2,3,4 




^ Figure 2.- Load distri- 

Figure Lift curve c> "bution of a 

* , °? + £ f a „ ta ~ tapered wing with twist. Figura 

pered wing with and ,ASUr,? 



without twist. 



3.- Wing with 
cutaway. 
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Fige. 6,6,7,8 
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Figs. 9,10,11,12,13 
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Figure 12.- Elliptical 

zero distri- 
bution and corresponding 
twist as used in the ex- 
ample. 




4*> 



0,50 



Flight condition 
c a total - +1.0 




Figure 13.- Results of lift distribution 
in first and fourth approach. 



.A.O.A. Teohnioal Memorandum No. 778 



Figs. 14,15,16,17,18 




Figure 14.- Elliptioal aero 

distribution in 
roll and separation into mean 
angle of roll and elliptical 
angle of twist in roll. 



-Theor. ourve, 
Blrnbaum-Glauert 
-Evaluation test 
data 




Figure 15.- Theoretical and 

experimental el- 
evator effect within range 
of small control deflections. 





Figure 16.- Angle of twist in 
roll and zero dis- 
tribution in roll of aileron. 




Figure 18.- Oomparison of data 
on normal distribu- 
tions of the Flavag III wing. 




Figure 17.- First and fourth approach 
of aileron calculation. 
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Figs. 19,20,21,22,33,34 




Figure 19.- Comparison 
of data on 
twisting effect for the 
Flavag III wing. 
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Figure 30.- Comparison 

of true * 
of both check calcula- 
tions. 
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Figure 31.- Comparison 

of true 
chord distribution based 
on results of lift dis- 
tribution analysis. 



to 



0.5 





















t 


j 


































X 














































'fix 


Y 


















7 


















/// 












































I 1° 





Qf 02 03 OA 05 



Figure 23.- Rolling mo- 
ment coeffi- 
cient of rectangular 
wing with aileron accord- 
ing to Lotz (ZW 1931). 
The small circles are those 
of Lippisch obtained as 
first approach. 
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Figure 23.- Comparative 
calculation 
of a zero distribution. 
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Figure 24.- Induced yawing 
moment coeffi- 
cient of rectangular wing 
with aileron (Lotz ZFM 1931) , 
The small circles are those 
of Lippisch obtained as 
first approaoh. 



